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Abstract 

The capacitated vehicle routing problem (CVRP) [TV02] involves distributing (identical) items 
from a depot to a set of demand locations, using a single capacitated vehicle. We study a general- 
ization of this problem to the setting of multiple vehicles having non-uniform speeds (that we call 
Heterogenous CVRP), and present a constant-factor approximation algorithm. 

The technical heart of our result lies in achieving a constant approximation to the following TSP 
variant (called Heterogenous TSP). Given a metric denoting distances between vertices, a depot r 
containing k vehicles having respective speeds {Ai}*L]^, the goal is to find a tour for each vehicle 
(starting and ending at r), so that every vertex is covered in some tour and the maximum completion 
time is minimized. This problem is precisely Heterogenous CVRP when vehicles are uncapacitated. 

The presence of non-uniform speeds introduces difficulties for employing standard tour-splitting 
techniques. In order to get a better understanding of this technique in our context, we appeal to ideas 
from the 2-approximation for scheduling in parallel machine of Lenstra et al. [LST90] . This motivates 
the introduction of a new approximate MST construction called Level-Prim, which is related to Light 
Approximate Shortest-path Trees [KRY95] . The last component of our algorithm involves partitioning 
the Level-Prim tree and matching the resulting parts to vehicles. This decomposition is more subtle 
than usual since now we need to enforce correlation between the size of the parts and their distances 
to the depot. 
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1 Introduction 



The capacitated vehicle routing problem (CVRP) is an extensively studied combinatorial optimization 
problem (see e.g., the book [TV02] and references therein). CVRP is defined on a metric space (V,d), 
where F is a finite set of locations/vertices and d : V x V ^ M_|_ a distance function that is symmetric 
and satisfies triangle inequality. There is a depot vertex r £ V that contains an infinite supply of an 
identical item, and each vertex u £ V demands some units qu of this item. A single vehicle of capacity 
Q > is used to distribute the items. The objective is to find a minimum length tour of the vehicle 
that satisfies all demands subject to the constraint that the vehicle carries at most Q units at any time. 

CVRP is closely related to the Traveling Salesman Problem (TSP). It is clear that CVRP reduces to TSP 
in the absence of capacity constraint. More interestingly, a reverse relation is also known — essentially 
the best known approximation algorithm for CVRP [HK85] achieves a guarantee of p+ 1, where p is the 
best approximation ratio for TSP. 

In practice, it is natural to have a fleet of multiple vehicles that can run in parallel. The objective can 
then be to either minimize the sum of completion times of all the vehicles or to minimize the maximum 
completion time over all vehicles (or the makespan of the routing). Furthermore the vehicles can all be 
identical (same speed) or heterogeneous (have different speeds). In either case, it is not hard to see that 
the total completion time objective reduces to the usual CVRP on a single maximum-speed vehicle, and 
constant-factor approximation algorithms readily follow. 

When the objective is to minimize the makespan with identical vehicles, ideas for approximating the 
regular CVRP problem using a tour-splitting heuristic introduced by Frederickson et al. [FHK78] can 
be easily adapted to derive a constant-factor approximation algorithm (see below). 

This motivates the Heterogenous Capacitated Vehicle Routing Problem (HetCVRP) that we consider. In 
this problem, a fleet of k vehicles with non-uniform speeds and uniform capacities is initially located at 
the depot vertex r. The objective is to satisfy the demands subject to the capacity constraints while 
minimizing the makespan. Our main result is a constant-factor approximation algorithm for HetCVRP. 

Most of our algorithmic ideas lie in solving the special case of HetCVRP when there is no capacity 
constraint. This problem, which we call HetTSP, is a generalization of TSP that might be of independent 
interest. For most of this paper, we will focus on obtaining a constant-factor approximation for HetTSP. 

1.1 Previous Techniques 

Tour-splitting solutions: To illustrate the use of known techniques, we outline how to obtain a 
constant-factor approximation algorithm for HetTSP with uniform speeds [FHK78]. First, notice that 
the union of the tours of OPT connects all vertices, and hence a minimum spanning tree has length at 
most k ■ OPT. Then consider an MST, duplicate its edge and take an Euler tour C, which is of length 
d{C) < 2k ■ OPT. Now split C into k segments of lengths at most by removing edges. Finally, 
the tour for the i^^ vehicle is obtained by connecting both endpoints of the i*'* segment of C to the 
depot. Since twice the distance from the depot to any vertex is a lower bound on OPT, the length of 
each tour is at most 3 • OPT and hence this solution is a 3-approximation. We remark that this can be 
extended to obtain an 0(l)-approximation for HetCVRP with uniform speeds (e.g., using Theorem 4.1 
in Section 4). 

At a very high level, this strategy has two main components: (1) Partitioning an MST into manageable- 
sized connected parts; (2) assigning these parts to vehicles. This simple idea — which was already 
present in the 70 's — is the central piece of many heuristics and approximations for vehicle routing 
problems (e.g., [FHK78, HK85, EGK+04, AHL06, GNRIO, AG87, AG90]). However, it is not clear how 
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to employ this technique in the presence of vehicles with multiple speeds. This is because the two main 
components now need some correlation: a small part of the MST, which should be assigned to a slower 
vehicle, must also be relatively closer to the depot in order to be reachable by this vehicle. 

Set-cover based solutions: For HetTSP with non-uniform speeds, previous approaches seem to give 
only a logarithmic approximation, as follows. Guess the optimal makespan OPT (within a constant 
factor). If each vehicle of speed s is given a length budget of s ■ OPT, then the vehicles can collectively 
cover all vertices. Using an approximation algorithm for A;- MST [Gar05] (or the related orienteering 
problem [BCK+07, CKP08]) within a maximum-coverage framework (see e.g.. [CK04]), we can obtain 
tours of length OPT that cover a constant fraction of all vertices. Repeating this coverage step until all 
vertices are covered gives a solution to HetTSP of makespan 0(log n) - OPT. The intrinsic problem of this 
approach is that it is too general — in fact, the above algorithm also yields a logarithmic approximation 
even in the setting where the metric faced by each vehicle is arbitrary (instead of just scaling by its 
speed), and this generalization of HetTSP can be shown to be set-cover hard. It is unclear whether the 
approximation of this set-covering based approach can be improved for HetTSP. 

1.2 Results, Techniques and Outline 

We extend the first tour-splitting approach described above to obtain the following result. 

Theorem 1.1 There are constant-factor approximation algorithms for HetTSP and HetCVRP. 

In order to obtain the approximation for HetTSP, we abstract the requirements of the two components 
in the tour-splitting strategy. As a preprocessing step, we round the speeds of vehicles to powers of two 
and guess the optimum makespan M. First, we specify conditions which guarantee that a collection 
of r-rooted trees is "assignable", that is, each vehicle can visit the nodes of the trees assigned to it 
within time 0{M) (Definition 3.1). The conditions in Definition 3.1 are based on the LP to obtain a 
2-approximation for scheduling in unrelated parallel machines by Lenstra et al. [LST90]. 

Secondly, instead of partitioning an MST as in the previous section, we consider more structured span- 
ning trees which we call Level-Prim trees. Consider grouping the vertices ordered according to their 
distance from r into levels, where the ith level includes all vertices within distance 2*M.^ The Level- 
Prim tree is simply the tree resulting from running Prim's algorithm with the restriction that all nodes 
in a level are spanned before starting to pull in nodes from the next. 

A Level-Prim tree has two important properties: (i) The vertices along every root-leaf path are mono- 
tonically nondecreasing in level and (ii) For every suffix of levels, the subgraph induced on it costs at 
most 0(1) times its induced MST. The first condition, which is the departing point from MSTs, greatly 
simplifies the decomposition procedure carried in the next step. The second property is related to the 
assignability conditions in Definition 3.1 and guarantees that the we can decompose a Level-Prim tree 
into an assignable collection. These properties are formalized in Theorem 3.3. 

The Level-Prim construction combine both MST and shortest-path distances from a root, so it is not 
surprising that this structure is related to Light Approximate Shortest-Path Trees (LAST) introduced by 
Khuller et al. [KRY95]. Indeed, we use the existence of a suitably defined LAST in proving Theorem 3.3. 
We remark, however, that the properties guaranteed by LASTs are not enough for our purposes (see 
Section 3.2). 

The third main component of our approximation for HetTSP is decomposing Level- Prim into an assignable 
collection of r-rooted trees. Roughly, we partition the edges of Level-Prim into subtrees while ensuring 

^Notice that given the rounding of vehicle speeds to powers of two, vertices in level i can only be served by vehicles of 
speed 2' or higher given the makespan bound M. 
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that each subtree consisting of vertices in levels up to i (and hence is at a distance of about 2*M from 
the root) also has length approximately 2*M, and thus can be assigned to a vehicle of speed about 
2*. This partition, which relies on the two properties of Level- Prim, gives a collection of unrooted trees 
which is assignable. Due to the length of these trees, the extra distance to connect them to the root r 
can be charged to their edges, hence this collection can be turned into a r-rooted assignable collection. 

In order to obtain an approximation to HetCVRP, we reduce this problem to approximating HetTSP in 
a suitably modified metric space. This new distance function encodes any additional trips to and from 
the root that a vehicle has to make if it runs out of capacity. The exact transformation is presented in 
Section 4. 

1.3 Related Work 

For the CVRP, the best known approximation ratio [HK85] is essentially p + 1 where p is the best 
guarantee for TSP. The current best values for p are p = \ for general metrics [Chr76], and p = 1 + e 
(for any constant e > 0) for constant dimensional Euclidean metrics [Aro98, Mit99]. This has been 
improved slightly to 1 + p • (1 - ^) - 3^ when Q > 3 [BDO07]. Recently, Das and Mathieu [DM10] 
gave a quasi-polynomial time approximation scheme for CVRP on the Euclidean plane. 

Several variants of TSP have been studied, most of which have a min-sum objective. One related problem 
with min-max objective is nurse station location [EGK"''04], where the goal is to obtain a collection of 
trees (each rooted at a distinct depot) such that all vertices are covered and the maximum tree length is 
minimized. Even et al. [EGK+04] gave a 4-approximation algorithm for this problem. This is based on 
partitioning the MST and assigning to trees along the lines of Section 1.1; their second step, however, 
involves a non-trivial bipartite matching subproblem. 

In proving the properties of Level-Prim, we use Light Approximate Shortest-Path Trees introduced by 
Khuller, Raghavachari and Young [KRY95], building on the work on shallow-light trees of Awerbuch, 
Baratz and Peleg [ABP90]. An (a,/3)-LAST is a rooted tree that has (a) length at most (3 times the 
MST and (b) the distance from any vertex to the root (along the tree) is at most a times the distance 
in the original metric. Khuller et al. [KRY95] showed that every metric has an ^a, 1 -|- ^^^^-LAST (for 
any a > 1) and this is best possible. 

One phase of our algorithm uses some ideas from scheduling on parallel machines [LST90], which also 
has a min-max objective. In this problem, job j has processing time pij on machine i and the goal is to 
assign jobs to machines while minimizing the maximum completion time. Lenstra et al. [LST90] gave 
an LP-based 2-approximation algorithm for this problem. 

2 Model and Preliminaries 

The input to the Heterogenous TSP (HetTSP) consists of a metric {V^d) denoting distances between 
vertices, a depot r ^ V and k vehicles with speeds {Aj}^^^ greater than or equal to 1. The goal is to 
find tours {rj}^^;^ (starting and ending at r) for each vehicle so that every vertex is covered in some 
tour and which minimize the maximum completion time max^^-^ 

At the loss of a factor of two in the approximation, we assume that the Aj's are all (non- negative integral) 
powers of 2. Then, for each integer z > we use pi to denote the number of vehicles with speed 2*. We 
let OPT denote the optimal value of this modified instance of HetTSP. 

We let G = {V, E) be the complete graph on vertices V with edge-weights corresponding to the distance 
function d. For any set F <Z E of edges, we set d{F) = Y2eeF^e- Given any (multi)graph H and a 
subset U of its vertices, H[U] denotes the subgraph induced on U and H/U denotes the graph obtained 
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by contracting vertices U to a single vertex (we retain parallel edges). Moreover, for any pair of vertices 
u, V in H, we use du^u, v) to denote the length of the shortest path in H between u and v. 

3 Algorithm for HetTSP 

Assume that we have correctly guessed a value M such that ^ < OPT < M. (This value can be found 
via binary search and we address this in the end of Section 3.3.) We partition the set of vertices V 
according to their distance to r: 

Vo = {ueV : d{r, u) < M}, and 

Vi = {ueV : d{r, u) G {2^~^M, TM]}, for aU i > 0. 

The vertices in Vi are referred to as level i vertices. For any i > 0, we use V<i as a shorthand for U*-^qV^- 
and similarly V<j = U*~qV^' = V<j_i. 

We also define the level of an edge {u, v) £ E as the larger of the levels of u and v. For each i > 0, Ei 
denotes the edges in E of level i. Note that de < 2^~^^M for all e G Ei, since both end-points of e are 
in V<i and the triangle inequality bounds its distance by the two- hop path via the root. We use the 
notation E<i = Uj^qEj and E>i = ^j>iEj. 

3.1 Assignable Trees 

We start by studying collections of trees that can be assigned to vehicles in a way that each vehicle 
takes time 0{M) to visit all of its assigned trees. 

Definition 3.1 (Assignable Trees) A collection of r -rooted trees IJi>o'^ covering all vertices V is 
called (a, /3) -assignable if it satisfies the following properties. 

1. For each i > and every T G 71, d{T) < a 2* M . 

2. For each i > 0, Y.j>i d(.Tj) < /3M Ej>i-i 2^ Mj- 

Intuitively, the trees in 71 can be assigned to vehicles with speed 2* so as to complete in time 0{aM). 
Condition (2) guarantees that the trees Uj>j '71 targeted by vehicles of speed 2* and above stand a 
chance of being handled by them within makespan 0{(5M). Interestingly, these minimal conditions are 
enough to eventually assign all trees in collection to vehicles while guaranteeing makespan 0{{a + I3)M). 

Lemma 3.2 Given an assignable collection IJi>o'^ of r -rooted trees, we can obtain in polynomial time 
an (4a -|- 2 P)- approximation for HetTSP. 

To prove this lemma^, we show that condition (2) guarantees the existence of a fractional assignment of 
trees where each vehicle incurs load at most /3M. Then using condition (1) and a result on scheduling 
on parallel machines [LST90], we round this assignment into an integral one while increasing the load 
on each vehicle by at most 2aM. We loose an extra factor of 2 to convert the trees into routes. 

■^We remark that a direct proof of Lemma 3.2 is also possible, but the route we take reveals more properties of the 
requirement at hand and could potentially be useful in tackling generalizations of HetTSP. 
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Fractional Assignment. Consider the bipartite grapli H whose left side contains one node for each 
tree in (J^ 71 and whose right side contains one node for each vehicle. (We identify the nodes with their 
respective trees/ vehicles.) There is an arc between the tree T G % and a vehicle of speed 2-^ if j > i — 1. 

Consider the following 6-matching problem in H: for each tree T S 71, we set b{T) = d(T) and for 
each vehicle u of speed 2^ we set b{u) = /32-^M. A (left-saturating) 6- matching is one which fractionally 
assigns all b{T) units of each tree T such that no vehicle u is assigned more than b{u) units. Notice that 
a feasible 6-matching gives a fractional assignment of trees where each vehicle incurs load at most /3M. 

Then our goal is to show the existence of a 6-matching in H. Using a standard generalization of Hall's 
Theorem (e.g., see page 54 of [CCPS98]), we see that H has a feasible 6- matching iff for every set V 
of trees, Y^tgv^C^) most J2ueN{v)^(''^)^ where N{V) is the neighborhood of V. However, the 

structure of H allows us to focus only on sets V which are equal to Uj>j 7j for some z.^ Using this 
revised condition, H has a fe-matching iff for all i, X^j>j d{Tj) < f3M X]j>j_i 2-'/Xj. Since this is exactly 
condition (2) in Definition 3.1, it follows that H indeed has a 6-matching (which can be obtained in 
polynomial time using any maximum flow algorithm [CCPS98]). 

Scheduling Parallel Machines. We show how to round the fractional assignment obtained in the 
previous section. We consider each tree as a "job" and each vehicle as a "machine" , where the "processing 
time" pt,u of a tree T in a vehicle u of speed 2^ is d{T)/2^; then the "makespan" of a vehicle is exactly 
equal to the sum of the processing times of the trees assigned to it. 

Let 

xt,u denote the fraction of tree T assigned to vehicle u given by scaling down a 6-matching in H 
(i.e., if the matching assigns d units of T to vehicle u, we have xt,u = d/d{T)). The feasibility of the 
matching gives YIt 

xt,uPt,u ^ (3M for all u. Moreover, by construction of the edges of H , xt^u > 
for r G 7i implies that u has speed at least 2*~^. Then using property (1) of assignable trees we get 
that xt,u > implies pT,u < 2aM. These two properties guarantee that x is a feasible solution for 
the natural LP formulation for the scheduling problem with a feasible makespan value of /3M and the 
maximum processing time t set to 2aM. Theorem 1 of [LST90] then asserts that x can be rounded into 
an integral assignment of trees to vehicles such that the load on any vehicle is at most (2a + /3)M. 

As in Section 1.1, we can transform each tree in Uj>o 71 into a cycle while at most doubling its length, 
which then gives a (4a + 2/3) approximation for HetTSP. This concludes the proof of Lemma 3.2. 

3.2 Level-Prim 

In order to obtain an assignable collection of r-rooted trees for our instance, we formally introduce 
Level- Prim trees. These are the trees obtained by the the following procedure. 

Algorithm 1 Level-Prim(G) 
1: For each i > 0, let Hi be an MST for G[V<i\/V<i. 
2: return V. = Ui>o^i- 



Note that Level-Prim trees can alternately be defined by modifying Prim's algorithm such that nodes 
in level i are only considered to be added to the tree after all nodes in levels below i have already been 
added. 

Theorem 3.3 A Level-Prim tree % = {-ffi}i>o satisfies the following: 

^To see that all other inequalities are dominated by those coming from such sets, first notice that if V contains a tree 
in Ti then NiV) already contains all vehicles of speed 2-' for j > i — 1. Then adding to V extra trees in Uj>i 7; does not 
change its neighborhood and thus leads to a dominating inequality. 
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• The vertex-levels along every root-leaf path are nan- decreasing. 

• For each i > 0, Y.j>i d{Hj) < 8 • MST {G/V<i). 

Note that the second property in Theorem 3.3 mirrors the second property in Definition 3.1. A formal 
connection between the two is estabhshed via the fohowing lemma that uses an optimal vehicle routing 
solution to derive a feasible spanning tree connecting a suffix of the level sets. 

Lemma 3.4 (Lower Bound) For each level £>0, MST(G/y<^) < M ■ J2j>e-i 2^ l^j- 

Proof: Consider an optimal solution for HetTSP and let E* be the set of edges traversed by vehicles 
in this solution; label each edge in E* by the vehicle that traversed it. Clearly E* connects all vertices 
to the root r. 

Observe that only vehicles having speed at least 2^~^ can even reach any vertex in V>£ (since a vehicle 
of speed s travels distance at most s ■ OPT < s ■ M). Thus every edge in E* n E>£ must be labeled by 
some vehicle of speed at least 2^~^. This implies that d {E* n E>£) < M • Ylj>e~i 2"' /^ji since the right 
hand side is a bound on the total length traversed by vehicles having speed at least 2^~^. 

On the other hand, since E* connects all vertices, E* D E>£ contains a spanning tree of G/V^e- Thus 
we have MSJ{G/V<e) <d{E*n E>i) < M ■ J2j>i-i 2^' ■ 
We then get the following corollary of Theorem 3.3. 

Corollary 3.5 A Level-Prim tree % = {-f^i}i>o satisfies the following: 

• The vertex-levels along every root-leaf path are non- decreasing. 

• For each i > 0, Y.j>i d{Hj) < 8M^j.>._i 2^ fij. 

In the rest of this section, we prove Theorem 3.3. It is easy to see that for every i, IJj=i spans 
G[y<£], hence the procedure produces a spanning tree for G. Moreover, by construction we obtain that 
every root-leaf path in T-L traverses the levels in non-decreasing order as desired. Thus, we focus on 
proving the second property in the theorem. 

Instead of comparing the length of the edges in T-L with an MST, it turns out to be much easier to use 
a specific LAST tree as proxy for the latter. The following LAST is implicit in the construction given 
in [KRY95]; for completeness we outline a proof in Appendix A. Recall that a spider is a tree with at 
most one vertex (the center) having degree greater than two. 

Theorem 3.6 ([KRY95]) Given any metric {V,d) with root r, there exists a spanning spider L with 
center r such that: 

• For each u £ V, the distance from r to u in L is at most 2 • d(r, u) . 

• The length of L is at most four times the MST in {V,d), i.e. d{L) < 4 • MST. 

We remark that we cannot use a LAST directly instead of Level-Prim since the former does not need to 
have the properties asserted by Theorem 3.3; it is easy to find a LAST which does not satisfy the first 
property, while Figure 1 also shows that the second can also be violated by an arbitrary amount. Using 
these spider LASTs we can obtain the main lemma needed to complete the proof of Theorem 3.3. 

Lemma 3.7 For any graph G and any Level-Prim tree % on G, we have d{'H) < 8 • MST(G'). 

Proof: Consider a spider LAST L for G and let V denote the set of all root-leaf paths in L; note that 
V is edge-disjoint. 
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(a) MST on G 



(b) A (1,2)-LAST T on G 



Figure 1: Instance with depot plus n nodes in Vq (box and gray nodes) and n nodes in Vi (black nodes). 
The distance between each distinct pair or nodes in Vq is 1, the distance from a gray node to a black 
node is also 1 and the distance between two distinct black nodes is 1/n. Notice that d(T D Ei) = n 
while MST(G/yo) < 2. 

Consider any root-leaf path P = (r = ui — )• M2 — • • • — ^ "Ufc) in "P. We claim that P crosses levels almost 
in an increasing order. Specifically, there does not exist a pair of nodes Ui, Uj G P with i < j, Ui £ Vi 
and Uj G V<e-2- Suppose (for a contradiction) that this were the case; then we would have that 

dL{r,Uj) = dL{r,Ui) + dL{ui,Uj) > dL{r,Ui) > d{r,Ui) > 2^~^M, 

where the last inequality uses Ui G Ve. On the other hand, d{r,Uj) < since uj G V<i-2\ so we 

obtain dL{r,Uj) > 2 • d{r,Uj), which contradicts the definition of L (see Theorem 3.6). 

Now we transform L into another spider L' which traverses levels in non-decreasing order as follows. 
For each root-leaf path P = {r = ui —7-^2 — s- ... Uf^.), perform the following modification. Let 
{ai,a2, ■ ■ ■ ,ak'} be the subsequence of P consisting of the vertices in even numbered levels, i.e. each 
G L2i for some £ > 0. Similarly, let {61, 62, • • • 1 ^fc"} be the subsequence of P consisting of the vertices 
in odd numbered levels. Define two paths Peven := (r —)• ai —)•... —)• a^i) (shortcutting P over nodes 6j's) 
and Podd = (r 61 bk") (shortcutting P over a^'s). Observe that both Peven and Podd cross 

levels monotonically: if not then there must be some i < j in P with Ui £ Vi and uj £ V<^_2, contrary to 
the previous claim. Also, by employing the triangle inequality we have that d{Peven), d{Podd) < d{P). 
Finally define the spider V as the union of the paths {Peven-, Podd} over all root-leaf paths P £ V. 

By construction, the vertex levels along each root-leaf path of L' are non-decreasing. Additionally 
d{L') = J2pevid{Peven) + d{Podd)) < 2 J2pevd{P) = 2 • d{L) < 8 • MST, by Theorem 3.6. Now 
partition the edges of L' as: 

r L'[Vo] if£ = 0, 

\ L'[V<e]\L'[V<e-i] i > I. 

By the monotone property of paths in L' , it follows that L'[V<e] is connected for every i > 0. Thus 
Ai is a spanning tree in graph G[y<^]/V<£. Since Hi in the Level-Prim construction, is chosen to be an 
MST in G[V<e\/V<i, we obtain d{Hi) < d{Ae). So, dCH) = Y.t>o d{He) < J2e>o = d{L') < 8 • MST. 
This completes the proof of the lemma. ■ 

Completing proof of Theorem 3.3. We now prove the second property in Theorem 3.3. Lemma 3.7 
directly implies this property for i = 0. For any level i > consider the graph G' = G/V^i] observe that 
\j-y^Hj is a Level-Prim tree for G' (due to the iterative construction of "H = |J^>q//^). Thus applying 
Lemma 3.7 to graph G' and its Level-Prim Uj>i Hj, we have Y.j>i d{Hj) < 8-MST("G') = 8-MST(G/V<i). 
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3.3 Decomposition Procedure 

In this section we decompose a Level- Prim tree into an assignable collection IJi>o^ r-rooted trees. 
Motivated by Corollary 3.5, the idea is to essentially break each subgraph Hi into many pieces and 
connect them to r in order to form the set of trees %. More specifically, assmne for now that each 
connected component in Hi is large enough, i.e. has length at least 2*M. Then for each i > 0, break 
the connected components of Hi into trees of length approximately 2'M; add to each tree the shortest 
edge connecting them to r and set 71 as the collection of r-rooted trees obtained. By construction we 
get that Ui>o satisfies the first property of an assignable collection. Moreover, notice that each edge 
added to connect a tree to the root has approximately the same length as the tree itself; this guarantees 
that d{7i) < 2d{Hi). It then follows that the collection Uj>o'^ assignable. 

Notice that it was crucial to break Hi into trees of length at least approximately 2*. But this is prob- 
lematic when Hi has a small connected component. In this case we show that such a small component 
is always attached to ( "dangling" from) a large enough component in (otherwise the dangling edge 
to a much earlier level will already make this component heavy enough not to be small); then we simply 
treat the small component as an integral part of the latter. 

Now we formally describe the proposed decomposition procedure. 

Step 1. Let Sq contain the subtree Hq = ^{^Eq. For each level i > 1: partition edges 'H{^Ei into 
a collection Si of (unrooted) subtrees such that each subtree contains exactly one edge from V<^i to Vi. 
For any r E 5j call the unique edge from V<j to Vi its head-edge h{T). Note that such a partition is 
indeed possible since '}i\y<i]/V<^i is connected. 

Any subtree in Si (for i > 0) is referred to as a level i subtree. Note that head-edges are defined only 
for subtrees in level 1 and above. 

Step 2. For each level i > 0: mark those t £ Si that have d{T) > 2^~^M. In addition, mark the tree 
Hq in Sq. Let S^ and 5" denote the marked and unmarked subtrees in Si. 

Step 3. For each level i > 1 and unmarked a G 5": define Tr{a) as the subtree in IJj^iSj containing 
the other end-point of h{a). 

Claim 3.8 For i > 1 and unmarked a G Sf , 7r(cj) G Moreover, TT{a) is marked. 

Proof: Since a is unmarked in level i > 1, d{h{a)) < d{a) < So the end-point v of h{a) 

in 7r(cr) satisfies d{r,v) > | • 2^~^M, otherwise d{h{a)) > 2*~^M - d{r,v) > 2*~^M. In particular 
V G Vi-i and thus vr(cr) G Si-i. 

For the second part of the claim, notice that if i = 1 then 7r(o") = Hq, which is always marked. 
So suppose i > 2. From the above, tt{<j) is in level i — 1 > 1 and hence contains a head-edge. 
This implies that tt{(t) contains some vertex u G y<i_i, namely an end-point of h('K{a)). But then 
d{TT{a)) > d{u,v) > d{r,v) — d{r,u) > 2*~^M, where we used d{r,u) < 2*~^ since u G V<i-i and 
d{r, u) > I • 2*~^M from above. Thus T^icr) must be marked. ■ 

Step 4. For each level i > and marked r G 5™: define Dangle(r) = 7r^-^(r) as the set of all unmarked 
a G ^.'Yi having 7r(cr) = r. Clearly d{a) < 2'-'^M for all a G Dangle(T). 

Step 5. For each level i > and marked r G Sl^: partition the tree r U Dangle(r) into subtrees 
Ti, . . . ,Tq such that the first q — 1 trees have length in the range [2'^~^^M, 2*+^M] and Tg has length at 
most 2*"'"^M. Notice that this is possible since all edges of rU Dangle(T) belong to i?<j+i and hence have 
length at most 2'^~^^M. Finally, add the shortest edge from r to each of these new subtrees to obtain a 
collection Ti{T) of r-rooted trees. 
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Claim 3.9 For any T e Tiir), we have d{T) < 3 • 2'+^M. 

Proof: Notice that every T G 7i(r) consists of a Tj (for some 1 < j < q) and an edge from r to 
a node in l/<j+i. Since the former has length at most 2*+^M and the latter has length at most 
2*+^M, it follows that d{T) < 3 • 2'+'^M. U 

Claim 3.10 Erer,(r) d{T) < 5 • [d{T) + (i(Dangle(r))]. 

Proof: We break the analysis into two cases depending of q. Suppose q = 1, namely Ti^r) 
consists of a single tree T. In this case T = r U Dangle(T) U {e}, where e is an edge to r. If 
1 = then d{e) = and the result holds directly. If i > then r has a node in V<^i and hence 
d(e) < 2'-^M. Because r is marked and different than Ho, the lower bound on its length implies 
that d{e) < 2'-^M < 4d(r) < 4(d(r) + d(Dangle(r)). The result follows by adding the length of 
T U Dangle(r) to both sides. 

Now suppose q > 1. Since all trees in 7i(r) lie in V<i+i, each edge from the root in Ti{T) has length at 
most 2*+^M. So the left hand side is at most T.'j=i d{Tj) + q ■ 2'+'^ M . But for j < q we have d{Tj) > 
2*+^M, so the last term of the previous expression can be upper bounded by ^^li-j Yl'jZi d{Tj). 
This bound is smallest when q = 2, which then gives ^xeTiir) '^(^) — •^Yl'j=i d{Tj) < 3d(T). This 
concludes the proof of the claim. ■ 

Step 6. For each level i > 0: define % = (Jre^™ Ti{T)- 

The following lemma summarizes the main property of our decomposition procedure. 

Lemma 3.11 The collection {Ti}i>o obtained from the above procedure is {Q, AO) -assignable. 

Proof: By Claim 3.9, each tree in % has length at most 3-2'''*"^M. So the collection satisfies condition (1) 
of Definition 3.1. 

Fix any i > for condition (2) in Definition 3.1. Due to Corollary 3.5, it suffices to prove that 
Yl,j>i d{Tj) < 5 • J2j>i d{Hj). Using Claim 3.10 we obtain that 

diTj) = diTiir)) < 5 • ^ [d{T) + d(Dangle(r))] = 5 • ) + 5 • d{SJ+,). 

The last equality above uses the fact that that {Dangle(T) : r G S-"} is a partition of Thus: 

J2 diTj) <5-Y^ d{Sf) + 5 • ^ d(5]Vi) < 5 • = 5-5] d{H,). (1) 

j>i j>j j>i j>j j>i 

This concludes the proof of Lemma 3.11. I 

Summary of the Algorithm. Our algorithm starts with an initial low guess of M and runs the Level- 
Prim procedure. If the second condition in Corollary 3.5 does not hold for this run, we double the guess 
for M and repeat until it is satisfied (this happens the first time that M reaches the condition for the 
correct guess: ^ < OPT < M). We use the decomposition in this section summarized in Lemma 3.11 
to obtain a (6,40)-assignable collection of trees. Using Lemma 3.2 on this collection gives us the desired 
constant approximation ratio by observing that the guess M in this step obeys M <2 ■ OPT. 
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4 Generalization for HetCVRP 



The input to the Heterogenous CVRP (HetCVRP) consists of a metric {V, d) denoting distances between 
vertices, depot r £ V (containing an infinite supply of items), demands {qv}v£V and k vehicles with 
speeds {Xi}^^i, each having capacity Q. A solution to HetCVRP consists of tours {ai}^^^ (starting and 
ending at r) for each vehicle so that all demands are satisfied and each vehicle carries at most Q items 
at any point in time. The objective is to minimize the maximum completion time, max^^^^ '^^y^ . We 
study the "split-delivery" version of CVRP here, where demand at a vertex may be served by multiple 
visits to it; however, our result easily extends to the "unsplit-delivery" HetCVRP. 

We show that the HetCVRP problem can be reduced to HetTSP in an approximation preserving way; so 
we also obtain an 0(l)-approximation for HetCVRP. The idea in this reduction is to modify the input 
metric based on lower-bounds for CVRP [HK85]. In order to avoid ambiguity, we use OPTyrp to denote 
the optimum for HetCVRP and OPTtsp to denote the optimum for HetTSP. 

Theorem 4.1 Consider an instance I of HetCVRP. There is a poly-time constructive instance J of 
HetTSP such that 0PTtsp(c7) = 0(1) • OPTvrp(X). Moreover, a solution to J of makespan M can be 
converted in poly-time to a solution to X with makespan 0{M). 

Proof: Let X be an instance of HetCVRP as specified above. Standard scaling arguments can be used 
to ensure that Q is polynomial in n and G {0, 1, . . . , Q} for aW. v (details in the full version). 

Let G = {y, E) denote the complete graph on vertices V with edge-weights equal to distances d. 
Augment G to a new graph H by adding vertices V = {vp : t> E y, p G [Qv]}-, and edges E' = {{v, Vp) : 
V £V,p £ [Qv]}] each edge {v,Vp) has weight '^^g^^ . For any vertex v &V, the vertices {vp : p G [iv]} 
are referred to as copies of v. Let iV' ,t) denote the metric induced on vertices V where d. denotes the 
shortest-path distances in graph H. We let J be the instance of HetTSP on metric {y\l) with depot 
r and k vehicles having speeds {Aj}^^]^. Since Q < poly{n) this reduction runs in polynomial time. 

For any graph L and subset S of vertices, let MinStL(S) denote the minimum length Steiner tree 
connecting 5. For any subset T C V and v €z V let Njj{T) denote the number of f -copies in T; 
also define it{T) = {v G V : N^iT) > 0}. Observe that for any T C V' we have MinStH(r) = 
MinStG(vr(r)) + Y.v&v ^"^^'''^^ by the definition of graph H. 

We first show that 0PTtsp(i7) = O(0PTvrp(X)). Consider an optimal solution {ai}^^^ to X. For each 
i G [k], let Ci{v) G [qv] denote the units of demand at vertex v G V served by vehicle i, and let 
Si = {v ^ V : Ci{v) > 0}. Note that Yli=i Ci(^^) = Qv for all v ^ V; hence we can choose S'^ C V' for each 
i G [k] such that U^^^S'^' = V and Ny{Sl) = Ci{v) for all w G i G [k]. Since cji is a capacitated tour in 

G serving demands {ci{v) : v G Si\, we have d{ai) > max |MinStG({r} U Si), YlveSt | using 

the (connectivity and capacitated routing) lower-bounds for CVRP [HK85]. Thus M'\nStH{{r} U S[) = 
MinStG({r} U Si) -\- Ylv&Si '^'''^^Q^'''^'' — ^ " '^i'^i)- consider the solution to J where the i*'* vehicle 

visits vertices S[ along the minimum TSP tour on {r} U 5^', for all i G [k\, the distance traversed by the 
z*^ vehicle is at most 2 • MinSt//({r} U 5"^') < 4 • d{ai). So the HetTSP objective value of this solution is 
at most maxjg[;j] = 4 • OPTyrp (2^). 

Now consider a solution {ri}^^^ to with makespan M. Let Ri C V' denote the vertices that are 
served by each vehicle i £ [k]. Since Tj is a TSP tour on {r} U Ri, we have d{Ti) > MinSt//({r} U Ri) = 
MinStG({r} U n^Ri)) + E^eV ^ ^ t^l and consider the instance of CVRP on 

vertices {r} U vr(i?j) with demands {N^{Ri) : v G 7r(i?j)}. As mentioned in the previous paragraph, 

max I l\/linStG({?'} U 7r(i?j)), J2veTr{Ri) ^"^'^q'^*'^'^^ i is a lower-bound for this instance, and the algorithm 
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from [HK85] returns a solution fij within a p = 0(1) factor of this lower-bound. It readily follows that 
{a}^^^ is a feasible solution to I with makespan at most maxjg[,t] = 0{M). ■ 

We note that this algorithm returns a non-preemptive HetCVRP solution, i.e., each item once picked up 
at the depot stays in its vehicle until delivered to its destination. Moreover, the lower-bounds used by 
the HetCVRP algorithm also hold for the (less restrictive) preemptive version, where items might be left 
temporarily at different vertices while being moved from the depot to their final destination. Thus our 
algorithm also bounds the "preemption gap" (ratio of optimal non-preemptive to preemptive solutions) 
in HetCVRP by a constant. 

5 Open Problems 

One interesting open question regards the approximability of HetTSP and HetCVRP when vehicles are 
located in multiple different depots across the space. The current definition of an assignable collection 
and the definition of Level-Prim crucially depend on the assumption of a unique depot, hence an extension 
to the multi-depot case is likely to require new ideas. Another interesting direction is to consider 
HetCVRP with non-uniform capacities, where the ideas presented in Section 4 do not seem to generalize 
directly. 
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A Proof of Theorem 3.6 

We will show that the following algorithm produces an (a, /3)-spider for G. 



Algorithm 2 (a, /3)-Spider 
1: Consider an MST for G and traverse it in preoder to obtain a path = S = {r = ui ^ U2 ^ 

...^Un) such that d{S^) < 2 • MST. 
2: for i from 1 to n do 
3: if ds{ui,Ui) > ad{ui,Ui) then 
4: add (ui, Ui) to S and mark Ui. 
5: for each marked node Uj, remove (nj_i,Mj) from 5. 
6: return S* = S. 



Lemma A.l The graph S* returned by the algorithm is a spider. 

Proof: Since the algorithm keeps adding edges to the path S^, it is clear that before Step 5 only the 
root ui and marked nodes have degree larger than 2. Moreover, each marked node has degree exactly 
3. Thus, after Step 5 we have that only the root has degree larger then 2, and the lemma follows. ■ 

Lemma A. 2 S* is an {a, (3) -spider. 

Proof: First we prove that ds*{ui,Ui) < ad{ui,Ui) for all i. To see this, consider S right before Step 
5. It follows from Step 4 that ds{ui,Ui) < ad{ui,Ui) for all i. Noticing that S* is a shortest path tree 
of S from node ui implies the desired result. 

Now we prove that S* satisfies the second property of (a, /3)-spider. Define vq = ui and let Vi be the ith 
node marked by the algorithm. It is clear that d{S*) < d{S^) + X]i=i d{ui,Vi); so our goal is to upper 
bound the last summation. 

Fix a node fj. Consider the beginning of the iteration where Vi is marked. Notice that at this point 
ds{ui,Vi) < d{ui,Vi^i) + ds{vi^i,Vi), since edge {ui,Vi^i) was already added to S; since is subgraph 
of S, it is also clear that the right hand side is at most d{ui,Vi^i) + dsQ{vi-i,Vi). However, since Vi 
was marked, we have that dsiui,Vi) > a ■ d{ui,Vi), and then using the previous bounds we obtain that 
a ■ d{ui,Vi) < d{ui,Vi-i) + dso{vi_i,Vi). 

Adding the previous inequality over all Uj's we get that a d{ui, Vi) < J2i d{ui,Vi-i) + '^^ dso{vi^i,Vi). 
Noticing that d{ui,vo) = d{ui, ui) = and reorganizing leads to (a — 1) d{ui,Vi) < dso{vi-i,Vi). 
Finally, notice that X^^Li dso{vi-i,Vi) < d{S^): this follows from traversing the path and using the 
triangle inequality. This gives the bound J2i d{ui, Vi) < d{S^)/{l — a). 

Plugging this back to a previous bound on the length of S* gives d{S*) < (1 + 1/(1 — a))d{S^) < 
(2 + 2/(1 — a))d{S^). This concludes the proof of Theorem 3.6. ■ 
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